Dependence of the intrinsic spin Hall effect on spin-orbit interaction character 
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We report on a comparative numerical study of the spin Hall conductivity in two-dimensions 
for three different spin-orbit interaction models; the standard k-linear Rashba model, the k-cubic 
Rashba model that describes two-dimensional hole systems, and a modified k-linear Rashba model 
in which the spin-orbit coupling strength is energy dependent. Numerical finite-size Kubo formula 
results indicate that the spin Hall conductivity of the k-linear Rashba model vanishes for frequency 
uj much smaller than the scattering rate r~^, with order one relative fiuctuations surviving out to 
large system sizes. For the k-cubic Rashba model case, the spin Hall conductivity does not depend 
noticeably on lot and is finite in the dc limit, in agreement with experiment. For the modified k- 
linear Rashba model the spin Hall conductivity is noticeably lot dependent but approaches a finite 
value in the dc limit. We discuss these results in the light of a spectral decomposition of the spin 
Hall conductivity and associated sum rules, and in relation to a proposed separation of the spin Hall 
conductivity into skew-scattering, intrinsic, and interband vertex correction contributions. 

PACS numbers: 72.10.-d,73.21.-b,73.50.Fq 
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I. INTRODUCTION 

Interest in spintronicsiiSiSi^i^ has been heightened by 
the technological impact of ferromagnetic metal based 
devices and by ferromagnetic semiconductor materials 
advances. Theoretical attention has recently focused on 
spintronics effects in paramagnetic materials, and in par- 
ticular on the spin Hall effeclii in which an electric field 
induces a transverse spin current. Murakami et al& and 
Sinova et al& have argued in different contexts that the 
spin Hall conductivity can be dominated by a contribu- 
tion that follows from the distortion of Bloch electrons 
by an electric field and therefore approaches an intrinsic 
value in the clean limit. The intrinsic spin Hall conduc- 
tivity adds to the skew scattering contribution that had 
been the focus of earlier theoretical work^LiSiii and can 
be altered by disorder vertex corrections. 

The proposed intrinsic spin Hall conductivity has 
drawn theoretical attention^^i"i^'^i^'^i^^i^^i^^i^^i^° to this 
unfamiliar transport coefficient. It has been argued 
that the intrinsic spin Hall conductivity does not sur- 
vive in the diffusive transport thermodynamic limit, 
either generally or for the specific case of the two- 
dimensional electron system with Rashba spin-orbit in- 



teractions (R2DES)*i studied by Sinova et 



Sev- 



eral researchers have pointed out that in the Rashba 
2DEG case, ladder vertex disorder correctionsi^siSiSS in 
the Kubo formula lead to a vanishing spin Hall cur- 
rent. The same conclusion has been reached via quantum 
Boltzmann theorjii^ii^ calculations which capture the 
same physics. These vertex correction claim is specific to 
the R2DES case, and has specifically been discounted for 
two-dimensional hole gasesSSiS^, p-doped bulk semicon- 
ductors and the modified Rashba coupling case.^'' Sev- 
eral numerical studies basing on the Landuer-Buttiker 
approach in mesoscopic systems with leads also have been 
performed, but lack of clear trends in the size dependence 



of the spin Hall conductances they evaluataSS^^SiSZ,, and 
the possibility of edge effects near the contacts, cannot 
be connected to the possibility of the spin Hall effect in 
the thermodynamic limit. 

In this article we approach these issues numerically 
by evaluating the finite-size Kubo formula^®'^^"^*' for the 
spin Hall conductivity for two-dimensional electrons with 
three different spin-orbit coupling models. The motiva- 
tion for the models we have chosen follows from the pars- 
ing of spin-Hall conductivity contributions illustrated 
schematically in Fig. Quite generally the charge and 
spin-current operators are diagonal in Bloch wavevector 
in any model with a spin-independent impurity potential, 
but have matrix elements that are both diagonal and off- 
diagonal in band index. (The models we study have only 
two-bands and it is usually convenient to consider them 
as representing the two spin states of a spin-1/2 parti- 
cle.) The presence of diagonal matrix elements means 
that current will not decay in a perfect crystal and the 
longitudinal conductivity (and Hall conductivity) is con- 
sequently limited by disorder scattering of Bloch states. 

We start from a linear response Kubo formula 
approach 
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(1) 



where 



j" is the a component spin or charge current in direc- 
tion, and H. is the volume of the system. When disorder 
is treated perturbatively, the standard ladder diagram 
approximation expresses the kernel K in Eq. (j^J in the 
following form.^'' : 
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J2 tr[j^G(k, C„ + u„,)J,G{k, Cn)l (3) 
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FIG. 1: Schematic diagrammatic representation of intrin- 
sic, interband coherence vertex correction, and skew scat- 
tering contributions to the spin Hall conductivity. The to- 
tal spin current (left vertex - open circle corresponding to 
= {v,Sz}/2) induced by an electric field (right vertex - 
dark circle corresponding to j = — ev) can be separated into 
a contribution from the density-matrix response that is off- 
diagonal in band index and a contribution from the density- 
matrix response that is diagonal in band index. The diago- 
nal response is dominated by a skew scattering contribution 
that is in general proportional to the Bloch electron scatter- 
ing time T. For the models studied here, the spin-current 
operator is purely off diagonal in band index and the skew 
scattering contribution is absent. The off-diagonal response 
has a leading contribution, the intrinsic spin Hall conductiv- 
ity, that is completely independent of disorder and is property 
of the host Bloch bands. Vertex corrections to the intrinsic 
spin-Hall conductivity are also independent of r in the limit 
of weak scattering, but depend on both band-structure and 
disorder potential. The solid lines in the figure represent Born 
approximation Greens functions. 



where G is a Matsubara formalism Born-approximation 
Greens function, is the current operator rcnormalized 
by ladder-diagram vertex corrections, and all four quanti- 
ties are matrices in certain bases. In this approximation 
the longitudinal conductivity ends up being dominated 
by terms that are diagonal in band index at each current 
vertex and end up being proportional to the Bloch state 
scattering time t. This contribution captures the Boltz- 
mann theory physics in which the current is due to field 
induced changes in the occupation probabilities of Bloch 
states. It turns out that the spin Hall conductivity does 
in general have a corresponding contribution, but only if 
scattering violates the principle of microscopic reversibil- 
ity, i.e. only if the a ^ b scattering rate is not equal 
to & — > a scattering rate where a and b are Bloch states. 
Indeed this property is violated when scattering ampli- 
tudes are evaluated beyond the Born approximation, as 
illustrated schematically in Fig. ^ and the spin Hall con- 
ductivity ends up being proportional to r x S", where the 
skewness 5' is a measure of the violation of microscopic 



reversibility.-— The terminology we use here is borrowed 
from the theory of the anomalous Hall effect in ferro- 
magnetic metals and semiconductors, which is strongly 
analogous to the spin Hall effect in paramagnetic metals 
and semiconductors. Wc identify the contribution to the 
spin Hall conductivity that is proportional to r and due 
to changes in Bloch state occupation probabilities, i.e. to 
a response that is diagonal in both wavevector and band 
indices, as the skew scattering contributioniiMi For the 
models we study here the spin-current operator is purely 
off diagonal in band index, as we explain explicitly be- 
low. It follows that the skew scattering contribution to the 
spin-Hall conductivity vanishes for the models we study. 

We define the intrinsic spin Hall conductivity as the 
interband spin Hall conductivity cr(w) of a disorder free 
system, which has a finite value in general in the a; — > 
dc limit. This quantity is a property of the band struc- 
ture of the perfect crystal, hence the term intrinsic. The 
dc limit of the interband conductivity can, however, be 
altered by disorder even in the limit of arbitrarily weak 
disorder, r — > oo as illustrated schematically in Fig. ^ 
In perturbation theory the factors of associated with 
Born approximation disorder scattering can be canceled 
by factors of r associated with products of two-Greens 
functions that have the same band index. The evalua- 
tion of the disorder correction to the interband spin-Hall 
conductivity requires the evaluation of a ladder sum. The 
renormalized current with the ladder correction is given 
by solving the vertex equation, 

J(k;z,z')=j(k) 

+ E ^-^^^^^^^G{k\ z)J(k'; z, z')G(k', z'). (4) 
k' 

Explicit evaluation of this ladder sum correction requires 
some approximations and can normally be accomplished 
only for very simple disorder models, or in the limit of 
small spin-orbit interactions. These limitations of per- 
turbation theory motivate the numerical study reported 
on here. 

We have reported previously on the influence of disor- 
der on the spin Hall conductivity of the k-linear Rashba 
model, concluding that it remains finite in the thermody- 
namic limit 1^ This conclusion is at odds with our current 
numerical findings, extrapolating to infinite system size 
and then to zero frequency, which are consistent with 
the perturbation theory conclusion that the dc spin Hall 
conductivity of this model is zero. The numerical stud- 
ies are complicated by the strong frequency dependence 
and large fluctuations in the spin Hall conductivity that 
occur in flnite-size calculations. The earlier calculations 
erred by using the frequency dependence of the longitu- 
dinal conductivity, which has corrections that vary like 
(tjr)^ compared to the (wr)^ dependence of the spin Hall 
conductivity discussed below, to judge whether or not 
the dc limit has been reached. The new findings super- 
sede the conclusions reached in Ref. (2^ with regards to 
the thermodynamic dc limit. Our numerical results for 
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the k-linear and k-cubic Rashba models are now consis- 
tent with analytic calculations which consider only the 
thermodynamic limit within the diffusive regime, hence 
resolving the controversy that has been associated with 
the linear Rashba model. '^^ 

In the k-cubic Rashba model, which approximately de- 
scribes two-dimensional valence band holes in a narrow 
quantum well with structural inversion asymmetrySS*^, 
the spin Hall conductivity does not show measurable ujt 
dependence. For this model, its a; — > limit is consistent 
with the pure intrinsic value. This observation suggests 
that the spin Hall induced edge spin accumulationsSSi 
recently seen in two-dimensional hole systems"^^ follows 
from the intrinsic spin Hall effect. We also study a modi- 
fied Rashba interaction which combines elements of the k- 
linear and k-cubic Rashba model and provides an approx- 
imate model for conduction band quantum well states in 
inverted gap materials like HgTe.'^^ As we explain below 
the modified model has spin-orbit splitting that varies 
like k^, as in the k-cubic Rashba model, but a wavevec- 
tor dependent Zeeman field whose in-plane orientation 
rotates once when the wavevector rotates once around 
the Fermi surface as in the k-linear Rashba model. In 
perturbation theory, the second property implies that an 
angular integral that appears in the vertex correction cal- 
culation and (for short-range impurity scattering) van- 
ishes in the k-cubic Rashba model, is non-zero. Vertex 
corrections to the intrinsic spin Hall conductivity survive 
for the modified model. Our numerical results demon- 
strate, however, that the corrections are present but do 
not cause the total spin Hall conductivity to vanish as it 
does for the k-linear Rashba model. The special prop- 
erty of the k-linear Rashba model that causes the spin 
Hall conductivity to vanish is related to the equation of 
motion of the spin-operatoriiiSS 

Our paper is organized as follows. In Section II we de- 
scribe the spin-orbit coupling and disorder terms in the 
model we study numerically. The spin-orbit interaction 
can be described in terms of a position and momentum 
dependent Zeeman field, whose orientation variation as 
a function of wavevector plays the key role in spin Hall 
conductivity calculations for these models. The disor- 
der model we employ assumes a scalar random poten- 
tial. We argue that as long as the random potential is 
dominantly spin- independent, this assumption is not es- 
sential. In Section HI we introduce the finite-size-system 
Kubo formula which expresses the spin Hall conductivity 
in terms of Hamiltonian eigenstates of a finite size two- 
dimensional electron system with area and periodic 
boundary conditions. The spin Hall conductivity eval- 
uated using this formula tends to fluctuate wildly from 
disorder realization to disorder realization and is very 
sensitive to avoided level crossings that occur close to 
the Fermi energy. These fluctuations are conveniently 
mitigated by evaluating the spin Hall conductivity for 
a continued complex frequency uj ^ z — irj along the 
imaginary axis. The dc spin Hall conductivity should be 
evaluated by first letting L'^ —> oo and then 77 —s- 0. Our 



expectation is that for systems much larger than a mean- 
free-path in size, dependence will appear only for 77 
smaller than or comparable to the finite size level spac- 
ing SE. Thus we should be able to extrapolate to the dc 
value as long as systems sizes can be reached numerically 
that are large enough to make other characteristic energy 
scales like the spin-orbit splitting and the life-time broad- 
ening energy h/r, much larger than rj 3> SE. In Sec. IV 
we present our numerical results for the finite-size Kubo 
formula and discuss its extrapolation to infinite system 
sizes. In Section V we discuss a spectral representation 
for the spin Hall effect and some associated sum rules. 
In Sec. VI we discuss the equation of motion of the spin 
operator for the k-linear Rashba model. Using the fact 
that the time derivative of the spin operator is propor- 
tional to the spin Hall current for this model, we are able 
to demonstrate that the finite size spin-Hall conductivity 
is always zero when averaged over boundary conditions. 
The typical size of the spin-Hall conductivity fluctuation 
in a given finite size system, is however, much larger than 
the intrinsic spin Hall conductivity. Finally in Sec. VII 
we present our conclusions. 



II. MODEL HAMILTONIAN 

We base our studies on a series of models with gener- 
alized Rashba spin-orbit interaction of the form: 

H = h^k'^/2m + iXg{\]i\){kZ<J+ - k°^a-) (5) 

where m is the effective carrier mass, k± — k^ i iky, 
and g is either unity or a function of |k|. (We set 
h = 1 for simplicity.) The Rashba model for an in- 
version asymmetric conduction band quantum well is 
generated by choosing a = I and g = 1. Valence 
band quantum wells have a more complex structured 
In the thin quantum well limit, light hole bands be- 
come energetically irrelevant and the heavy hole bands 
can be effectively described by Eq.(l) with a = 3, the 
so called k-cubic Rashba model. We study in addition 
a modified Rashba model with g — k^ and a = 1 (i.e. 
Hb — Afc^[k X z] • o"), which provides an approximate 
model for conduction electron quantum well states in in- 
verted gap materialsi^ In the following we use the Fermi 
energy Ep and the inverse Fermi wave number k^^ in the 
absence of both spin-orbit coupling and disorder as the 
units of energy and length respectively. We take a dis- 
order model consisting of uncorrelated short-range scalar 
impurity potentials: Viy) = ^^(r — R/) which sat- 

isfies {V{q)V{q^)*) = 7Viy2^(q- q'), where ^(q) is the 
Fourier component of V{y) and Ni — UiL^ is the num- 
ber of impurity scatterers which we take to correspond 
roughly to the number of carriers. We choose this type 
of disorder potential model, rather than the more real- 
istic finite correlation length model utilized in previous 
studies^^ in order to connect more directly with the an- 
alytical results. 
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We diagonalize the finite-size disordered electron 
Hamiltonian in the A = eigenstate basis and introduce 
a hard cutofi^ at a sufficiently large momentum A. (This 
means, of course, that the disorder potential has an ef- 
fective correlation length ~ A~^.) Our calculations are 
performed at a fixed carrier density Ue ~ kp /2Tr and at 
finite system sizes (see below) up to 70k p^, larger than 
the mean free path ^ lOkp^ and the Fermi wavelength 
^ kp . 

Finally we introduce the charge and spin current opera- 
tors. The charge current definition, j = ±ev = —dH/dA 
follows from the charge conservation continuity equation, 
where the vector potential A must be included to ob- 
tain a gauge invariant expression. In a system with 
a spin-rotational Hamiltonian charge spin components 
along arbitrary quantization axes are conserved sepa- 
rately and we can introduce a spin-dependent vector po- 
tential, A = (±e)Ac + Sq • Aq, the spin current is given 
by = —dH/dAa- With this definition a continuity 
equation expresses local conservation of each Cartesian 
component of spin. We retain the same definition when 
spin-orbit interactions are included, although the conti- 
nuity equation is now violated because the Hamiltonian 
is spin-dependent. 



Numerical evaluation of the spin Hall conductivity is 
complicated by the substantial fiuctuations in finite size 
system values when rj is small. Following the seminal 
arguments of Thouless and KirkpatrickS^, the physically 
appropriate value for 77 is 77 ~ gSE where g = 2EpT is 
the Thouless dimensionless conductance. The values of t 
quoted in our results were calculated from the golden-rule 
expression for the transport scattering rate, 

;i/r = 2^^|y(k-k')P(l-k-k')<5(^k' - ^f), (7) 

k 

which determines the Drude longitudinal charge conduc- 
tivity via (Jd — ne^r/m = 2EpT{e'^ /h). The variance 
of numerical spin Hall conductivities does appear to get 
smaller with system size, to the extent that this trend can 
be judged from our numerical results, but relative fiuc- 
tuations in magnitude arc still larger than one at small 77 
even for the largest system sizes that we are able to study. 
In our calculations, the disorder averaged spin Hall con- 
ductivity always has the same sign as the intrinsic spin 
Hall conductivity a^y, negative for the k- linear Rashba 
model, and positive for the k-cubic Rashba model. In 
Fig.(l) we have chosen a sign convention in which the 
intrinsic spin-Hall conductivity is defined as positive. 



III. FINITE SIZE KUBO FORMULA FOR 
SPIN-HALL CONDUCTIVITY 

We start from linear response theory, Eq.(l) and 
Eq.(2). An elementary calculation leads to the following 
formally exact expression for the static z-spin component 
spin Hall conductivity, 

^ »ft fjEg) - fjEg') {a\j^\a'){a'\jy\a) , , 

a, a' 

where = edH/dpi, = {dH/dpy, s^}/2, with being 
{h/2)az for electrons and {ih/2)crz for holes. In Eq. ® irj 
can be regarded as a complex frequency continued from 
the real axis to the imaginary axis and can be interpreted 
as an electric-field turn on time. In metallic systems, like 
the ones considered here, 77 must exceed the simulation 
cell level spacing 5E in order to obtain bulk values of the 
transport coefficients considered. At the same time, 77 
must be smaller than all other intensive energy scales 
such as the Fermi energy Ep, the spin-orbit coupling 
splitting A50, and the disorder broadening ft/r, where 
T is the scattering time. The finite value of 77 represents 
the coupling of a finite subsystem of a macroscopic con- 
ducting sample to its environment, leading for metallic 
systems to the loss of resolution of the discrete individual 
energy levels of the subsystem. For a finite system with 
periodic boundary conditions, the spin-Hall conductiv- 
ity is a function of 5E/Ep,r]/Ep,/^so/Ep, and H/rEp. 
The macroscopic dc spin Hall conductivity is obtained by 
extrapolating finite size results first to SE ^ (L — > 00) 
and then to 77 — > 0. 



IV. SPIN HALL CONDUCTIVITY NUMERICAL 
RESULTS 

A. Rashba model 

In Fig|21we plot spin-Hall conductivities evaluated for 
models with Epr/H = 8 with Xkp/Ep = 0.2 and various 
system sizes as a function of rj in the k-linear Rashba 
model (top panel), the k-cubic Rashba model (middle 
panel) and the modified Rashba model (bottom panel). 
We note that in the regime where 77 is smaller than 
Ep, Aso a-iid h/r but larger than SE (plotted as a func- 
tion of system size as an inset in Fig|2Jl, the spin Hall con- 
ductivity sometimes changes as a function of 7/, strongly 
so in the k-linear Rashba model case. Taking the limit 
77 — > 0, extrapolating from the regime where 77 > 6E 
is satisfied, we find that for the linear Rashba model 
the spin Hall conductivity is strongly suppressed in the 
thermodynamic limit. Our numerical results are consis- 
tent with the conclusion from analytic calculations that 
asH vanishes for this model. The disorder strength de- 
pendence of the spin-Hall conductivity for the k-linear 
Rashba model is shown in Fig. 3 where asH is plotted as 
a function of rjT and rj/Ep (inset) fixing rEp at 10, 8, 6, 
and 4. These results are consistent with the analytic the- 
ory conclusion that asH approaches the intrinsic value for 
Lu > r~^, but that it vanishes for this model for lu —>■ 0. In 
perturbation theory, the contribution to asH that varies 
on the frequency scale comes from vertex corrections 
to the intrinsic interband response. Our previous numer- 
ical results which reached an incorrect conclusion on the 
dc value of asH for this model, were performed at a value 
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FIG. 2: Spin Hall conductivity for a variety of system sizes 
as a function of rj in the k-linear Rashba model (top panel), 
k-cubic Rashba model (middle panel) and modified Rashba 
model (bottom panel). All of these results are for the case 
Eft/H = 8 and Xkp/EF ~ 0.2 Inset in panel a): Level spac- 
ing as a function of system size. The results shown in this fig- 
ure correspond to level spacings varying between ~ 0.003-Ef 
and ~ O.OOS-Bf. Inset in panel b): Disorder strength depen- 
dence of the spin Hall conductivity for the Rashba model. 
The dc (Limjj^oLimiB^ocsff (??, ^-E)) spin Hall conductivity 
appears to vanish for the k-linear Rashba model case only. 



of rj which gives accurate values for the dc longitudinal 
conductivity but, as we have now learned by extrapolat- 
ing ?7 ^ 0, not for the spin Hall conductivity. 



B. k"^ Rashba model 

For the Rashba 2DHG the spin Hall conductivity is 
insensitive to 77; there is no evidence of a relatively con- 
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FIG. 3: Spin Hall conductivity as a function of the relaxation 
time T in the k-linear Rashba model plotted as a function 
of r]T. Inset: Same quantity plotted as function of disorder 
strength Eft — 10, 8, 6, 4. These results demonstrate that 
the spin Hall conductivity has a contribution with frequency 
dependence on the scale r~^, and are consistent with the con- 
clusion tha the dc spin Hall conductivity vanishes. 



tribution that has frequency dependence on the scale 
of . This finding indicates that vertex corrections 
to the intrinsic spin Hall conductivity, at least for the 
short-range disorder scattering model we have studied, 
are absent or not notably large. We find similar behavior 
over a wide parameter range of A and r. This finding is 
consistent with analytical studies of the spin Hall con- 
ductivity by Bernevig and Zhang'^^ who find the vertex 
corrections vanish in this model. Our results should be 
compared with those from other recent numerical stud- 
ies of mesoscopic spin transport based on lattice versions 
of the Luttinger model for p-doped semiconductors. In 
both quantum well^ and bulk casesSS., the intrinsic spin 
Hall conductance is found to be robust against disorder 
in agreement with the present result. Also, as shown in 
the inset, the exact-diagonalization numerical results are 
in very good agreement with analytic calculations which 
correct for disorder only by including finite-lifetime cor- 
rections to the Greens functions that appear in the in- 
trinsic diagram. (This approximation is referred to in the 
figure as the relaxation time approximation.) This im- 
plies that the intrinsic effect will be dominant in systems 
with strong spin-orbit coupling satisfying Xkp ^ Ht^^. 



C. Modified Rashba model 

The modified Rashba model has properties interme- 
diate between those of the Rashba and k'^ Rashba 
models since the energy spectrum is identical to that of 
the k'^ model, but the eigenstates are the same as those 
of the k^ model. In this case the rj dependence of the 
spin Hall conductivity is shown in Fig. 12 (c). We find 
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a smooth deviation of cTsh at small 77 region. The nu- 
merical calculations do find indications of dependence on 
Tyr, implying that vertex corrections to the intrinsic spin 
Hall conductivity do not vanish. However, the overall 
spin Hall conductivity remains finite in the limit 77 — )■ 0. 

No analytic studies of vertex corrections in this model 
have been reported to date, although related models^ 
have been studied. Using a formalism similar to that 
developed in Ref. , we evaluate the vertex correction 
to the renormalized current. In the limit of A^qT ^ 1 
the vertex correction can be expressed in terms of 5ji, = 
Ji, — ji, = Sjlai. For short range scatterers. 



1 



^tr[a,G(k,z)j,G(k,^')], (8) 



where vp is the density of states in the absence of spin- 
orbit coupling. After a straitforward calculation, we ob- 



tain the following expression 5jl, 



vp^kp^] ~ {3/2h'pX)e^i[iyF-f-kpj^ + h'p-kp_]. where kp± 
and i'F± are the Fermi wave length and the density of 
state in the ± band. The spin-dependent Fermi wave- 
lenghts and densities-of-states that appear above are 
given by 



vf± = j^f(1 i SmXkp) 



(9) 



and 



kp- = — - 



2mA 



(1 - ^/l-SimXkp)^). (10) 



In the weak spin-orbit interaction limit, Ago/Ep << 1, 
these may be approximated by i>p± — > j^_f(1 =F SmXkp) 
and kp-i- ~ kp- —2'mXkp, and then vertex correction 
ends up being Sj = —Xkp[z x cr] which ends up canceling 
the intrinsic contribution to the spin-Hall conductivity. 
On the other hand, beyond the small spin-orbit coupling 
limit the exact cancellation between intrinsic and vertex 
contributions does not take place, consistent with the 
above numerical result. This is in sharp contrast with 
the k-linear Rashba model where the cancellation appears 
to hold for arbitrarily strong spin-orbit coupling. We 
comment further on this special property of the k-linear 
Rashba model later in the paper. 



V. SPECTRAL REPRESENTATION OF THE 
SHE AND SUM RULES 

The finite-size Kubo formula for the spin-Hall conduc- 
tivity, Eq.® may be expressed in the form 

h Im[(Q!|j,^|a')(a'|7x|a)l 



a, a' 



This form is based partly on our finding that the dissi- 
pative contribution to the spin Hall conductivity, which 
is not included in the above expression is vanishing, in 



analogous with the case of a charge Hall conductivity, the 
dissipative term strictly vanishes when spatial invariance 
is recovered by averaging over disorder realizations. 

It is instructive to consider the following spectral de- 
composition of the spin-Hall conductivity, 



dE 



NjE) 

E"^ + 77^ 



(12) 



where 



^(^) = jjil^[f{Eo)'f{Ea.')]lTnMfy\a'){a'\j^\a)] 

a, a' 



5{E~E^ + E^,). 



(13) 



In the following we first focus on the ordinary k-linear 
Rashba model. Ush depends on both the phase and the 
magnitude of the matrix elements in Eq.® and on the 
energy differences of the levels involved. The size of the 
matrix elements is characterized by the integral of N{E) 
over all energies which satisfies the following sum rule: 



dEN{E) = -^^/(i?„)Im[(a|[j-,j,]|a)] 



-hx 



a 

2mL2 



Y,!{E^){a\^y<7-\a) 



(14) 



The final form for the zeroth moment of the Hall spectral 
function {M'^) in Ea. (|14|l follows from the observation 
that the two terms in the Rashba spin-orbit interactions 
must have identical expectation values if isotropy is re- 
covered in the thermodynamic limit. For h/r Aso, 
(Hso) is close to its value in the perfect crystal state. We 
note that this expression is valid both in the presence (see 
below) and absence of an external magnetic field and that 
TTy is the kinetic momentum in the y direction. 

The left panel in Fig. 4 shows the spectral function 
N{E) as a function of E at Xkp/Ep — 0.2 and at two 
disorder strength EpT — 10 and EpT — 4. N(E) has 
a positive peak at E corresponding with the spin-orbit 
splitting 2Xkp {kp+ and kp- are approximately equal) 
and becomes negative with very small magnitude at small 
E. The negative contribution at small E corresponds to 
the vertex correction contribution to the spin Hall con- 
ductivity. The large peak near the spin orbit splitting 
energy corresponds to the intrinsic contribution to the 
spin Hall effect. In the limit of small spin-orbit coupling, 
the energetic width of the intcrband peak vanishes and 
the intrinsic spin-Hall conductivity is proportional to the 
ratio of the sum rule and the square of the spin-orbit 
splitting. The vertex correction contribution is enhanced 
in asH by small energy denominators. 

As first emphasized by Rashba«i2i, insight into the spin 
Hall conductivity of the linear Rashba model can be 
achieved by introducing an external magnetic field. Let- 
ting 



fik -ihV + eA(r), 



(15) 
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FIG. 4: The spin Hall spectral function N{E) in the absence 
(top left panel) and presence (bottom left panel) of a magnetic 
field. The right panels show the corresponding single-electron 
energy spectra in the absence of disorder. 



we introduce a magnetic field B perpendicular to the 
plane. (A(r) = BySc.) The charge current and the spin 
current in this case are given by j = —e{i/h)[H,r] = 
— e(7r/m — Xz x a) and = {■K/ni,az/2}/2, where 
TT = —iKV + eA(r) is the kinetic momentum. We con- 
sider the case of Landau level filling factor v ~ 1 and 
\/IbOJc = 2 as an example in the following. The en- 
ergy spectra are linear in the Landau level index in both 
bands and the spin-orbit splitting near the Fermi level 
is approximately bhujc as shown in right bottom panel 
in Fig. 4. We find that the spectral function N{E) has 
several peaks. The two left most peaks can be identified 
as the intra-band contributions that evolve into the ver- 
tex correction at zero field while the rest correspond to 
the inter-band contributions that give the intrinsic spin 
Hall effect. Interestingly the intraband contribution has 
both positive and negative peaks, with the negative peak 
appearing at lower energy and therefore having a larger 
contribution to the spin-Hall conductivity. As pointed 
out by Rashba^^ the intraband and interaband contribu- 
tions to asH cancel, as in the zero field case with disorder 
discussed above. 



VI. WHAT IS SPECIAL ABOUT THE 
K-LINEAR RASHBA MODEL? 

Perturbation theory calculations and these numeri- 
cal calculations consistently indicate that the k-linear 
Rashba model has zero spin Hall conductivity for any 
scalar random potential. This property is special to the 
k-linear model and is not generic, as indicated for exam- 
ple by our numerical results for the other models con- 



sidered in this paper. How should we understand this 
exceptional behavior? One possibility is to explain it in 
terms of vertex corrections to the model's current oper- 
ator, which cancel spin-dependent contributions on the 
Fermi surface. This explanation is not fully satisfactory, 
however, since the spin-Hall conductivity depends partly 
on contributions away from the Fermi surface, and the 
same cancellation does not cause the charge Hall con- 
ductivity of corresponding ferromagnet models to van- 
ish. Instead, as discussed previously elsewhere in the 
literaturei2iS£, the source of the special behavior is al- 
most certainly the following exact relationship between 
the time dependence of the total spin operator and the 
spin-current operator: 



•ru 1 /fe 2mA 



(16) 



(This type of relationship between in-plane spin equation 
of motion and in-plane spin-current, can be generalized to 
any model with spin-orbit coupling that is linear in mo- 
mentum, so that most of the conclusions reached below 
to the 2D Dresselhaus model and even to models that 
include both Dresselhaus and Rashba interactions and 
do not have circular Fermi lines.) Since the spin-density 
must approach a constant in the steady state, it follows 
that the non-equilibrium steady state expectation value 
of the left-hand-side of this equation, and hence also of 
the right hand side, must vanish. We believe that this 
argument is essentially valid, although there is subtlety 
in its application because the conductivity is defined by 
taking the thermodynamic limit and then the dc limit. 
(The same argument, naively applied, could be used to 
prove that the drift velocity vanishes in the steady-state 
induced by an external electric field.) Below we discuss 
the implications of this identity for the finite-size spin 
Hall conductivity calculations, and then discuss some 
of its implications for semiclassical descriptions of spin 
transport in the Rashba model and other systems. 



A. Finite-size Kubo formula 

We start by considering the linear response of an in- 
plane spin component to a constant change in the vector 
potential in a finite-size system. The vector potential 
change gives rise to a perturbation H' = j ■ 6 A: 



SSf, 



L2 ^ E, 

n.n' 



dt([s^(t),j,])e*('^+"')VL2 

{n\Sf_,\n'){n'\j^\n) 



fi Jo 

1 V- f{En)-f{E'J 



(17) 



For a square finite size system with side L, a change in 
vector potential by SA^, corresponds to a change in the 
boundary condition phase by 2TrLSA,^/^Q where $o is 
the electron flux quantum. Using the relation between 
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the time dependence of the spin operator and the spin- 
current operator, it follows that 

{i/h){En - E,,,){n\Sf,\n') = ~2m\{n\f^\n') /H^ . (18) 

Comparing with Eq. ^ we find that 



ti^ 6S, 



2mA SA 



(19) 



where z — uj + irj is a, complex frequency. This is an exact 
expression for the k-linear Rashba model with arbitrary 
scalar impurities. For a disordered R2DES s^; in general 
has a boundary condition and disorder potential depen- 
dent expectation value oc L, corresponding therefore to 
a spin-density per unit volume that vanishes in the ther- 
modynamic limit. Because the z — > value of the spin 
density response to vector potential is proportional to the 
derivative of the ground state spin density with respect 
to boundary condition phase angle in the direction of the 
vector potential, 



hm asH 

77^0 



2m\L 



(20) 



it is evident that the average of the z ^ value of this re- 
sponse function over boundary conditions is zero. (Since 
the spin-density must be a periodic function of (j)y with 
period 27r, the integral of its derivative over any period 
must vanish. In Ea. (|20|l (f)y is the boundary condition 
phase angle in the j/-direction and Sx = Tr[p^].) This 
appears to be the conclusion that can most confidently 
be drawn about Kubo formula properties from Eq. (|16f) . 
For T] — 0, the typical value of the spin Hall conductivity 
at a particular boundary condition is large in magnitude, 
indicative of large persistent spin currents in finite size 
systems. Our numerical results for the spin-Hall conduc- 
tivity appear to be consistent with the natural ansatz 
that averaging over boundary conditions is equivalent to 
averaging over disorder realizations in finite-system cal- 
culations of the spin-Hall conductivity. The fact that 
these averages at 77 = appear to yield the same val- 
ues for the spin-Hall conductivity as extrapolations from 
r/ > SE, guarantees that the equation of motion argu- 
ment for vanishing spin Hall conductivity in the linear 
Rashba model is valid. 

We note that it is possible to establish that the spin- 
Hall conductivity vanishes at integer Landau level fill- 
ing factors in the absence of disorder without appealing 
to a concrete calculation by using total spin equation of 
motion identities. In this case, the boundary condition 
phase angle 0j, just corresponds to an a;-direction guid- 
ing center shift in a translationally invariant system. In 
sharp contrast to the zero field case, there is no 0^ de- 
pendence in the absence of disorder; the result for any 
boundary condition equals the zero result obtained by 
averaging over boundary conditions. Because of the gap 
between Landau levels, no subtleties arise in taking the 
thermodynamic limit. The spin Hall conductivity clearly 
vanishes at any integer filling factor. 



There literature contains some arguments that the spin 
Hall conductivity vanishes for any model. For example, 
Sheng et alj^i have performed a numerical simulation 
similar with present work, and have concluded that the 
spin Hall effect of source of the zero spin Hall conduc- 
tivity of the linear Rashba model is more general. In 
particular they argue that because all energy eigenvalues 
have anti-crossing behavior as a function of boundary 
conditions (or equivalently flux <I> through a cylindrically 
shaped sample) spin transports cannot arise. We note 
that the linear response regime attains with external field 
eEL — —e d^/dt small compared with all relevant en- 
ergy scales but larger than the level spacing SE to gener- 
ate Landau-Zener tunneling through anti-crossing gaps. 
Consequently the adiabatic argument of Ref. 30] cannot 
capture the linear response of the spin transport. We 
rather conclude that strong suppression of the spin Hall 
conductivity is an accidental property of particular mod- 
els, not a generic effect. As we have seen in this article 
the k-cubic model gives a good example of a model for 
which the intrinsic spin Hall effect is dominant. The 
spin Hall effect observed in this model can be under- 
stood as an intrinsic eSeciM^. We note that the spin 
Hall effect is observed as spin accumulation. Although 
there is no analytic theory of spin accumulation due to 
an intrinsic spin Hall effect, we note that numerical stud- 
ies of spin accumulation^^'^^ show little accumulation in 
the k-linear Rashba model and robust accumulation in 
the k-cubic model. These results are consistent with the 
present spin transport study and a naive theory of spin 
accumulation. 



B. Implication for Semiclassical Theory of Coupled 
Charge and Spin Transport 

From the equation of motion for the density matrix 
it follows quite generally that the time derivative of the 
spin-density can be related to the equation of motion for 
the averaged spin-density: 



dt 



dt 



-Tr 



(21) 



In the steady state which balances acceleration by an ex- 
ternal electric field with disorder scattering, the density 
matrix p is constant and the spin-density is expected to 
saturate at a finite value. As discussed above, it follows 
from this argument that oc (d/(ii)Tr[pcr^/2] vanishes. 
It is interesting to compare the single-band semiclassi- 
cal theory of spin transport^^ with this result. In a ho- 
mogeneous system this theory describes the spin-density 
dynamics by the following equation: 



dS^ 
dt 



1 

Z2 



5]tr 
k 

2mA 



Q 

jz,int _ ^fjJ_ 



(22) 
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The first term in the first form for the right hand side cor- 
responds to the spin torque term which describes spin- 
density dynamics in the absence of colhsions. The ef- 
fect of collisions which scatter electrons between Bloch 
states is accounted for by the second term.^^ In the final 
form for the right hand side we have introduced the re- 
laxation time approximation for the scattering term and 
recognized that the collision free expression for the spin- 
density evolution is proportional to the Hall spin-current 
in the absence of collisions and hence to the intrinsic 
spin Hall conductivity. In the strong spin-orbit scatter- 
ing limit on which we focus, the spin relaxation timaii^ 
that appears in this equation may be approximately iden- 
tified with the momentum relaxation time. The steady 
state in-plane spin-density induced by an electric field is 
proportional to the intrinsic spin Hall conductivity and 
to the momentum relation time r. We know that above 
semiclassical argument of the spin-Hall conductivity fails 
to capture vertex corrections, because it does not prop- 
erly account for the influence of disorder on the interband 
components of the density matrix response^^ Since ver- 
tex corrections also change the value of the in-plane spin 
induced by an electric fieldji^ this theory may also fail 
to account quantitatively for the value of in-plane spin- 
density induced by an external electric field. 

VII. SUMMARY 

In this article, we have studied spin transport driven 
by an electrical potential bias in two-dimensional electron 
and hole systems with spin-orbit coupling due to struc- 
tural inversion asymmetry, primarily using finite-size ex- 
act diagonalization as a tool. We have studied three dif- 
ferent models of spin orbit coupling, the standard k-linear 
Rashba model, a k-cubic Rashba model appropriate for 
two-dimensional hole systems in narrow quantum wells, 
and a modified k-linear Rashba model with eigenspinors 
like that of the standard Rashba model and eigenval- 
ues like that of the k-cubic Rashba model. In these sys- 
tems a current of spins oriented perpendicular to the two- 
dimensional layer flows perpendicular to the direction of 
the electric field, an effect known as the spin Hall effect. 
For the models we have studied the expectation value of 
the spin current in the perfect crystal Bloch eigenstates 
is zero, implying that there is no skew-scattering-induced 
Bloch state occupation number change contribution to 
the spin Hall effect. The spin Hall effect is entirely due to 
interband coherence induced in the system by the elec- 
tric field. When disorder is treated perturbatively the 
spin Hall effect can be separated into an intrinsic contri- 
bution that is a property of the perfect crystal electronic 
structure along, and a disorder-related vertex correction 
contribution that remains finite even when the scattering 
rate vanishes. This vertex-correction is partially analo- 



gous to the scattering angle weighting correction that 
vertex corrections introduce into the theory of the longi- 
tudinal conductivity. 

The three models we study differ qualitatively on how 
vertex corrections alter the intrinsic spin Hall effect. We 
evaulate the spin Hall conductivity numerically for a fi- 
nite system and at a finite frequency it] continued to the 
imaginary axis. The frequency r/ can be thought of as 
a turn on rate for an electric field, or as energy level 
broadening due to the coupling of the small system be- 
ing studied numerically to the rest of a macroscopic sam- 
ple. The thermodynamic limit dc spin Hall conductivity 
must be calculated by first letting the system size be- 
come large and then letting 77 — > 0; the vertex correction 
appears as a dependence of the spin Hall conductivity 
on TjT. In the k-linear Rashba model we find that the 
spin Hall conductivity depends strongly on rjT, vanish- 
ing for ?7 — > 0. This finding is consistent with analyt- 
ical calculations have shown that the vertex correction 
strongly suppresses the intrinsic contribution to the spin 
Hall conductivity for this model. For the k-cubic Rashba 
model we find that vertex corrections vanish, a finding 
that may hold only for the short-range disorder scatter- 
ing model we apply. For the modified k-linear Rashba 
model, the vertex corrections do not vanish and alter 
the intrinsic spin Hall effect by a fraction that decreases 
with increasing spin-orbit coupling strength. Taken to- 
gether, these results demonstrate that the intrinsic in- 
terband spin Hall conductivity can be altered by vertex 
corrections, depending on details of the electronic struc- 
ture and the disorder potential. The special situation 
which leads to a vanishing total spin Hall conductivity for 
the k-linear Rashba model is related to the relationship 
between spin equations-of-motion and spin-currents that 
applies only for systems with spin-orbit coupling that is 
linear in momentum. Since the spin Hall conductivity of 
the k-cubic model is purely interband, and vertex cor- 
rections are weak for this model, we conclude that the 
spin-Hall-induced spin accumulation observed in a two- 
dimensional hole gas by Wunderlich et air- (for which 
the k-cubic model is applicable), must be due primarily 
to the intrinsic spin Hall effect. 



Acknowledgement 

The authors thank D. Culcer, T. Jungwirth, M. 
Koshino, S. Murakami, N. Nagaosa, Q. Niu, S. Onoda, 
E.I. Rashba, D. N. Sheng for useful discussions. One 
of the authors K.N. is supported by the Japan Society 
for the Promotion of Science by a Research Fellowship 
for Young Scientists. This work has been supported by 
the Welch Foundation and by the Department of Energy 
under grant DE-FG03-02ER45958. 



10 



^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman.J. M. 

Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka- 

nova, and D. M. Treger, Science 294, 1488 (2001). 
2 S. Parkin et al. Proc. of IEEE 91, 661 (2003). 
^ Semiconductor Spintronics and Quantum Computation, 

edited by D. D. Awschalom, D. Loss, and N. Sarmarth 

(Springer- Verlag, Berlin, 2002). 

* I. Zutic, J. Fabian, S. Das Sarma, Rev. Mod. Phys. 76, 
323-410 (2004). 

^ A.H. MacDonald, N. Samarth, and P. Schiffer, Nature Ma- 
terials 4, 195 (2005). 

" T. Dietl and H. Ohno, MRS BuUiten 28, 714 (2003). 
M.I. Dyakonov and V. I. Perel, Zh. Eksp. Ter. Fiz. 13, 657 
(1971) [JETP 33, 467 (1971)]. 

* S. Murakami, N. Nagaosa, and S. C. Zhang, Science 301, 
1348 (2003). 

® J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 
and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 
^° J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 
" S. Zhang, Phys. Rev. Lett. 85, 393 (2000). 

D. Culcer et al., Phys. Rev. Lett. 93 046602 (2004). 

J. Inoue, G. E. W. Bauer and L.W. Molenkamp, Phys. 
Rev. B 67, 033104 (2003); i6id 041303 (2004). 
" E.G. Mishchenko, A.V. Shytov, B.I. Halperin, Phys. Rev. 
Lett. 93, 226602 (2004). 
A. Khaetskh, cond-mat/0408136 

R. Raimondi, P. Schwab, Phys. Rev. B 71 033311 (2005). 
O.V. Dimitrova, cond-mat/0405339 

E. I. Rashba, Phys. Rev. B 70, 161201 (2004). 
" E.I. Rashba, Phys. Rev. B 70, 201309 (2004). 
^° O. Chalaev and D. Loss, cond-mat/0407342 

21 E. I. Rashba, Sov. Phys. Solid State 2, 1109 (1960). 

22 J. Schliemann and D. Loss, Phys. Rev. B 71, 085308 
(2005). 

23 B.A. Bernevig et al, Phys. Re v. B 70, 113301 (2004); B.A. 
Bernevig and S-C. Zhang, |cond-mat/0411457| 



2'' S. Murakami, Phys. Rev. B 69, 241202 (2004). 

2"^ L. Sheng, D. N. Sheng, and C. S. Ting, Phys. Rev. Lett. 



94 016602 (2005). 
2** B. K. Nikolic, L. 



Zarbo, and S. Souma, 



cond-mat/0408693 
2^ E. M. Hankiewicz, L.W. Molenkamp, T. Jungwirth, and J. 

Sinova, Phys. Rev. B 70, 241301 (2004). 
2* Thouless and Kirkpatric, J. Phys. C 14, 235 (1981). 
2^ K. Nomura, J. Sinova, T. Jungwirth, Q. Niu, and A. H. 

MacDonald, Phys. Rev. B 71 041304 (2005). 
3° D. N. Sheng, L. Sheng, Z. Y. Weng, F. D. M. Haldane, 

cond-mat/0504218 

31 K. Nomura et al., in preparation. 

32 B. K. Nikolic, S. Souma, L. P. Zarbo, J. Sinova, 
cond-mat/0412595 

33 A recent reexamination of the results of Mischenko et al 
14 by N. Sugimoto, S. Onoda, S. Murakami, N. Nagaosa, 
cond-mat/0503475 has challenged some of the conclusions 
reached in the quantum kinetic equation literature on the 
R2DES spin Hall effect. 

3" R. Winkler, Phys. Rev. B 62, 4245 (2000). 

35 Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 

Awschalom, Science bf 306, 1910 (2004). (2004). 
3^ J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 

cond-mat/0410295 
3'^ X.C. Zhang et al. PRB 63 245305 (2001). 
3* M.W. Wu and J. Zhou, cond-mat/0503616] 
3^ W.Q. Chen et al., cont-mat/0502570 

*° G. Mahan, Many-Particle Physics, 3rd edition, Kluwer, 

New York, 2000. 
*i P. Leroux-Hugon and A. Ghazali, J. Phys. C: Solid State 

Phys, 5, 1072 (1972). 
*2 N.A. Sinitsyn, Q. Niu, J. Sinova, K. Nomura, 

cond-mat/0502426 accepted for publication in Phys. Rev. 

B. 



